Anomalous Rashba spin splitting in two-dimensional hole systems 
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It has long been assumed that the inversion asymmetry-induced Rashba spin splitting in two- 
dimensional (2D) systems at zero magnetic field is proportional to the electric field that characterizes 
the inversion asymmetry of the confining potential. Here we demonstrate, both theoretically and 
experimentally, that 2D heavy hole systems in accumulation layer-like single heterostructures show 
the opposite behavior, i.e., a decreasing, but nonzero electric field results in an increasing Rashba 
coefficient. 
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Spin degeneracy of electron and hole states in a solid 
stems from the inversion symmetry in space and time. 
If the spatial inversion symmetry is broken, then there 
is a splitting of the single particle states even at mag- 
netic field B = (Ref. [I]). In quasi two-dimensional 
(2D) semiconductor structures, the bulk inversion asym- 
metry (BIA) of the underlying crystal structure (e.g., a 
zinc blende structure), and the structure inversion asym- 
metry (SIA) of the cppfining potential contribute to the 
B = spin splittingE While BIA is fixed, the so-called 
Rashba spin splittings due to SIA can be tuned by means 
of external gates that change the electric field E in the 
sample .DO The B = spin splitting is of significant cur- 
rent interest both because of its fundamental importance 
as well as its possible device applications.^ 

For many years it has been assumed that the Rashba 
spin splitting in 2D systems is proportional to the elec- 
tric field that characterizes the inversion asymmetry 
of the confining potential.Hu In single heterostructures, 
where SIA is the dominant source of spin splitting, the 
electric field is determined by the density-dependent 
self-consistent potential.!] One would thus expect that 
the spin splitting decreases with density, although for 
2D electron systems this effect may partly be compen- 
sated by many-particle effects that tend to increase the 
Rashba spin splitting for low densities.!!! Here we demon- 
strate, both theoretically and experimentally, that 2D 
heavy hole systems in accumulation layer-like single het- 
erostructures show the opposite behavior, namely, a de- 
creasing, but nonzero electric field results in an increas- 
ing Rashba coefficient. Contrary to electrons, however, 
exchange-correlation effects in the low-density regime de- 
crease the spin splitting. We show that this surprising 
result is a consequence of heavy hole-light hole (HH-LH) 
coupling in 2D hole systems. We obtain good qualitative 
agreement between calculated and measured spin split- 
tings in 2D hole systems in GaAs heterostructures where 
the density and the spin splitting are varied by means of 
an external gate. Our results are applicable to many sys- 
tems, as most III-V semiconductors have essentially the 
same band structure that is underlying our investigation. 

To lowest order of the wave vector k and electric field 



E, the SIA spin splitting of electron statesjji the T§ con- 
duction band is given by the Rashba ternxl 

fff c ° = ak x E-ct. (1) 

Here er = (a x , a y , a z ) denotes the Pauli-, spin matri- 
ces and a is a material-specific prefactor.El We assume 
E = (0,0, E z ). Treating the off-diagonal k • p coupling 
between electron and hole states by third order Lowdin 
perturbation theoryjj we obtain for the Rashba coefficient 
ct\ of the lowest electron subband A = 1 
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where P is Kane's momentum matrix elementS and 
A%(, = £1 - £$ with £1, £\, £[, and £{ the en- 
ergy of the Ath electron, HH, LH, and split-off subband, 
respectively.Eij The numerical prefactor a depends on the 
geometry of the confining quantum well (QW). In an in- 
finitely deep rectangular QW we have a = 256/(8l7r 2 ). 
According to Eq. (Q) we obtain a spin splitting ±aE z k\\ 
of the subband dispersion £ (k||) that is proportional to 
the electric field E z and is linear in the in-planc wave 
vector k|| = (k x ,k y ,0). A detailed analysis reveals that 
spin splitting of electron states depends on the electric 
field E v in the valence band that differs from the elec- 
tric field E c in the,|eonduction band by the contributions 
of the interfaces £3113 However, the important point here 
is that in a single heterostructure both E v and E c are 
determined by the self-consistent Hartree potential. 

For hole systems in the Tg valence band (point group 
Td), the dominant contribution to Rashba spin splitting 
is given by the terixo 

fff°=/?kxE-J, (3) 

where (3 is a system-dependent prefactor and J = 
{Jx, Jy, Jz) are the angular momentum matrices for j — 
3/2. We neglect here the small corrections in H^® due 
to the k • p coupling to remote-bands such as the higher 
Tg and Tj conduction bands Quantum confinement 
reduces the symmetry from the cubic point group T<z to 
D 2 d- Therefore, the four-fold degeneracy of the Tg band 
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(Td) is lifted and we obtain two- fold degenerate subspaces 
transforming according to T7 and Tq of the point group 
Did (Refs. |l6],|l7|) . The irreducible representation T-? cor- 
responds to the LH states (z component of angular mo- 
mentum j z — ±1/2) whereas Tq corresponds to the HH 
states (j z — ±3/2). To lowest order in k the effective 
Rashba Hamiltonian for the LH states is the same as in 
Eq. (||) for electron statesJla For HH states, however, spin 
splitting is mediated by a coupling to the LH states so 
that, to lowest order, spin splitting of HH states is of 
third order in k (Ref . |l4|) . Neglecting anisotropic correc- 
tions in the Hamiltonian we have 
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with a± = l/2(a x ± i(J v ) and k± = k x ± ik y . Treating 
the off-diagonal HH-LH coupling by third order Lowdin 
perturbation theorycl we obtain for the Rashba coefficient 
(3% of the lowest HH subband A = 1 
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well approximation!] we have, for the subband energies 
measured from the corresponding bulk band edge, 

El (x E* /3 which implies 8 V X cx A 2 / 3 and fl\ oc N~ 4 / 3 . 
Therefore, we can expect from Eqs. (||) and (||) that ac- 
cumulation layer-like 2D HH systems show a Rashba spin 
splitting that increases when N and E z are reduced. On 
the other hand, the coefficient (|2|) is essentially indepen- 
dent of N and E z because the energy gam,A^, are al- 
ways of the order of the fundamental gapO 

In order to validate these qualitative arguments we 
present next the results, of realistic, fully self-consistent 
subband calculations.EJ We use an 8 x 8 multiband 
Hamiltoniar£3 that includes the lowest conduction band 
Tg, the topmost valence band Tg and the split-off valence 
band T7. The simpler 4x4 Luttinger Hamiltonian,E£l 
taking into account only the band Tg, gives essentially 
the same results. We have checked that higher conduc- 
tion bands have a minor influence. Many-particle ef- 
fects are takea into account based on a density-functional 
approach. c3c3 From these calculations we obtain the dif- 
ference A A = N + — A_ between the spin subband densi- 
(5) ties N± as a function of the total density N = N + + N-. 
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where 72 and 73 are the Luttinger parameters .0 In an in- 
finitely deep rectangular QW we have a = 64/(97r 2 ). We 
see from Eq. (||) that the Rashba spin splitting of HH 
states depends not only on the electric field E z but also 
on the separation between the HH and LH subbands. A 
decreasing separation gives rise to an increasing Rashba 
coefficient The factor 73(72 ±73) in Eq. (||) refers to 
a quantum structure grown in the crystallographic direc- 
tion [001]. The expressions for other growth directions 
are similar, but the other terms in Eq. (|^) remain un- 
changed. We remark that for typical hole densities only 
the lowest HH subband is occupied. 

The electric field E z that enters into the Rashba Hamil- 
tonian depends on the charges in the system. We will 
show now that accumulation layer-like single heterostruc- 
tures behave rather differently with respect to changes of 
the 2D charge density as compared to other quasi 2D 
semiconductor structures. 

In a rectangular QW, a small density iV and a small 
asymmetry imply that the properties of the system are 
controlled by the effective potential steps at the inter- 
faces, i.e., changes in N or E z have a minor effect in 
this regime. In an inversion layer-like heterostructure, 
we always have a band bending of the order of the fun- 
damental gap so that, for small densities, the Hartree 
potential and E z are determined by the space charges 
due to the given concentration of ionized majority impu- 
rities in the system. For accumulation layer-like systems, 
on the other hand, it was shown by Sternal that the 
space charge layer is controlled by the much smaller con- 
centration of minority impurities in the system. Thus, 
even for a small 2D density, the dominant contribution 
to the Hartree potential stems from the charges in the 
2D system itself. Therefore, over a wide range of den- 
sities N, the electric field E z is proportional to N. In 
single heterostructures, the subband separations are ap- 
proximately proportional to E z . Using the triangular 
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FIG. 1. (a) Spin splitting AN/N and (b) effective spin 
splitting coefficient {PiE z )/{fj,h) as a function of iV for a 2D 
HH system in the accumulation layer of a GaAs-Alo.5Gao.5As 
single heterostructure on a (001) GaAs substrate, calculated 
including exchange-correlation (EXC, solid lines) and neglect- 
ing exchange-correlation (dashed lines). For the dotted lines 
see text. 

In Fig. |l|(a) we show AN/N calculated as a func- 
tion of N for a 2D HH system in the accumulation 
layer of a GaAs-Alo.sGae-sAs single heterostructure on 



a (001) GaAs substrateo From N — 5 x 10 11 cm~ 2 to 
1 x 10 10 cm~ 2 the parameter r 8 , the Coulomb energy 
to Fermi energy ratio, increases from 4.3 to 17. There- 
fore, one can expect that many-particle effects are quite 
important in this regime of densities N. Indeed, we 
find that taking into account exchange-correlation (solid 
lines) reduces AN/N as opposed to a calculation with- 
out exchange-correlation (dashed lines). This behavior, 
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which is opposite to 2D electron systems,@ can be traced 
back to the fact that exchange-correlation increases the 
subband spacingsE3 so that the Rashba coefficient f5\ is 
reduced, in agreement with Eq. (|B|). 

It is convenient to characterize our numerical results in 



terms of an effective Rashba coefficient (P^E Z ) (Ref. p5| ). 
Assuming that the spin-split HH subband dispersion is 
approximately of the form £±(k\\) = ± (J3iEg)k$, 

where fih (times 2/fi 2 ) is the reciprocal effective mass, we 
have 



($E Z ) = d-0i h ) 



N (^N + AN - - AN) + AN (y/N + AN + y/N - AN) 
6N 2 + 2AN 2 



(6) 



Figure |l|(b) shows that (/3{ l i? 2 )/(/i/ ! .) increases when N is 
reduced. 

For comparison, we have calculated AN/N for 
a 2D electron system in the accumulation layer 
a Gao.47Ino.53As-Alo.47Ino.53As single heterostructure 
[Fig. §(a)]. Here spin splitting is given by Eq. (@). There- 
fore, the spin-split subband dispersion is approximately 



of the form £±(k\\) — (Mc)&jj ± 
similarly to Eq. (0) (Ref. pq) 



aiE z )k^, and we obtain 



(aiE z ) = V2t7 <jt ) (VN + AN-VN- AN) . (7) 



In Fig. g(b) it can be seen that, in contrast to Fig. 
|f](b), the spin splitting coefficient (aiE z ) / (fi c ) decreases 
rapidly with decreasing N. We remark that unlike the 
HH system in Fig. |l|, exchange correlation has only a 
weak influence on the electron system in Fig. || (Ref. 27). 
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FIG. 2. (a) Spin splitting AN/N and (b) effective 
spin splitting coefficient (aiE z ) / (fj. c ) as a function of TV 
for a 2D electron system in the accumulation layer of 
a Gao.47Ino.53As-Alo.47Ino.53As single heterostructure (solid 
lines). For the dotted lines see text. 

To further analyze our numerical results, we define an 
effective electric field (E z ) = {d z Vn{z)), where Vh(z) is 
the Hartree potential without the effective potential due 
to the position-dependent band edges. In an accumula- 
tion layer the contribution of the space charge layer to 
Vh(z) is very small.E3 It follows then, by partial inte- 
gration of the Poisson equation, that (E z ) = e/ (2seo) N, 
where e is the dielectric constant and the expectation 



value refers to the 2D charge density that gives rise to 
Vh{z). Using these values for (E z ) and (fj, c ) = 89 eVA 2 
we obtain (cti) » 34.3 eA 2 independent of N, consistent 
with Eq. (g). This implies that in Fig. p|(b) the drastic 
change of (a%E z ) j (/i c ) merely reflects the change of the 
electric field (E z ). On the other hand, the weak varia- 
tion of (ffiE z )/{(j, h ) in Fig. [j](b) indicates that the "bare" 
Rashba coefficient (/3( l ) increases by a factor of 250 when 
N is lowered from 5 x I0 11 to I x I0 10 cm~ 2 (Ref. |§). 
This is in good, qualitative agreement with the analyt- 
ical model discussed above that predicts an increase of 
P'l by a factor of 50 4/3 w 184. Note that for low densi- 
ties the third order perturbation approach, that underlies 
Eqs. (^) and (|), breaks down because the subbands are 
merging together so that higher order corrections become 
important. These higher order terms are fully taken into 
account in our numerical calculations ffl If the density 
N is reduced below 10 10 cm~ 2 the Hartree potential and 
spin splitting are ultimately controlled by the fixed con- 
centration of minority impurities. c3c3 

It is interesting to compare the spin splittings in ac- 
cumulation layers with those in QWi's where E z is tuned 
externally, e.g., by means of gates.cl The dotted lines in 
Figs. |] and || show the calculated results for a 200 A wide 
rectangular QW where the external electric field E z xt was 
chosen according to E z xt (N) = e/(2eeo)N. In an elec- 
tron system (Fig. |^) we obtain spin splittings very close 
to the results for the accumulation layer. In particular, 
we have (011) ~ 30.6 eA 2 independent of N. Similarly, for 
a 2D HH system in a QW (Fig. 0) and N < 1 x f 11 cm" 2 
we obtain (P^) w 7.54 x 10 6 eA 4 . [For larger densities 
higher order corrections in £±(kn) become important Jiff] 
Since in QWs the subband spacings are essentially de- 
termined by the QW width (i.e., are independent of N) 
this is consistent with Eq. (H). These calculations also 
indicate that for 2D HH systems in a QW, spin spfLU 
ting becomes negligible in the regime of low densities ,E_3 
which is due to the fact that spin splitting of £±(k\\) is 
proportional to However, for 2D HH systems in sin- 
gle heterostructures, spin splitting can be very important 
in the low-density regime. We note that inversion layers 
give results similar to QW's, but the specific numbers 
depend on the details of the doping profile. 

In order to reinforce our conclusions, we present next 
a comparison between measured and calculated spin 
splittings in a GaAs-Alo.3Gao.7As single heterostructure 
grown on a nominally undoped (3ff)A GaAs substrate 
with a weak p-type background doping.Efl A back gate 
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was used to tune the density N from 1.8 x 10 to 



4.2 x 10 iu cm" 



To measure the spin subband den- 



sities N±, the Shubnikov-de Haas (Sd±D oscillations at 
low magnetic fields B were examinedotj at a tempera- 
ture T ~ 50 mK (see inset of Fig. ||). The frequencies 
/sdH of theses-oscillations are a measure of the zero-i? 
spin splitting^ In Fig. ||(a) we present the measured 
and calculated spin subband densities exhibiting remark- 
ably close agreement. Fig. 0(b) shows (0^ E z )/(nh) de- 
termined by means of Eq. (g). On average, (P^E z }/{/j,h) 
increases as the density is reduced. Taking into ac- 
count the orders-of-magnitude change that we have for 
E z ) I (fj,h) in QW's and for {a\E z ) /{fi c ) in electron 
systems, the agreement between experiment and theory 
is quite satisfactory.E3 We wish to emphasize that it is 
indeed the anomalous enhancement of the Rashba coeffi- 
cient in 2D HH systems in accumulation layer-like single 
heterostructures that allows us to experimentally resolve 
the spin splitting in this density regime; data on QW 
samples with comparable densities reveal no measurable 
spin splitting.Eil 
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FIG. 3. Measured (circles) and calculated (solid lines) spin 
subband densities N± (a) and effective spin splitting coeffi- 
cient ($E,)/(ji h ) (b) as a function of density N = N+ + N- 
for a 2D HH system at a GaAs-Alo.3Gao.7As single het- 
erostructure on a nominally undoped (311)A GaAs substrate 
with a weak p-type background doping. The inset shows the 
measured magnetoresistance R xx as a function of magnetic 
field B (upper part) and the Fourier transform (FT) of R xx 
(lower part) for N = 2.75 x 10 10 cm -2 . 
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